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The time-dependent Ginzburg-Landau approach is used to investigate nonlinear response of a 
strongly type-II superconductor. The dissipation takes a form of the flux flow which is quantitatively 
studied beyond linear response. Thermal fluctuations, represented by the Langevin white noise, are 
assumed to be strong enough to melt the Abrikosov vortex lattice created by the magnetic field into 
a moving vortex liquid and marginalize the effects of the vortex pinning by inhomogeneities. The 
layered structure of the superconductor is accounted for by means of the Lawrence-Doniach model. 
The nonlinear interaction term in dynamics is treated within Gaussian approximation and we go 
beyond the often used lowest Landau level approximation to treat arbitrary magnetic fields. The 
I-V curve is calculated for arbitrary temperature and the results are compared to experimental data 
on high-r c superconductor YBa2Cu307-«5. 
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I. INTRODUCTION 

Electric response of a high temperature superconduc- 
tor (HTSC) under magnetic field has been a subject of 
extensive experimental and theoretical investigation for 
years. Magnetic field in these layered strongly type-II 
superconductors create magnetic vortices, which, if not 
pinned by inhomogeneities, move and let the electric field 
to penetrate the mixed state. The dynamic properties of 
fluxons appearing in the bulk of a sample are strongly 
affected by the combined effect of thermal fluctuations, 
anisotropy (dimensionality) and the flux pinning 1 . Ther- 
mal fluctuations in these materials are far from negli- 
gible and in particular are responsible for existence of 
the first-order vortex lattice melting transition separating 
two thermodynamically distinct phases, the vortex solid 
and the vortex liquid. Magnetic field and reduced dimen- 
sionality due to pronounced layered structure (especially 
in materials like Bi2Sr2CaCuOs+5) further enhance the 
effect of thermal fluctuations on the mesoscopic scale. On 
the other hand the role of pinning in high-T c materials 
is reduced significantly compared to the low tempera- 
ture one, leading to smaller critical currents. At elevated 
temperatures the thermal depinning- further diminishes 
effects of disorder. 

Linear response to electric field in the mixed state 
of these superconductors has been thoroughly explored 
experimentally and theoretically over the last three 
decades. These experiments were performed at very 
small voltages in order to avoid effects of nonlinear- 
ity. Deviation from linearity however are interest- 
ing in their own right. These effects have also been 
studied in low-T c superconductors experimentally 2 ^ and 
theoretically^ and recently experiments were extended 
to HTSC compounds^. 

Since thermal fluctuations in the low-T c materials are 
negligible compared to the inter-vortex interactions, the 
moving vortex matter is expected to preserve a regu- 



lar lattice structure (for weak enough disorder). On the 
other hand, as mentioned above, the vortex lattice melts 
in HTSC over large portions of their phase diagram, so 
the moving vortex matter in the region of vortex liquid 
can be better described as an irregular flowing vortex 
liquid. In particular the nonlinear effects will also be 
strongly influenced by the thermal fluctuations. 

A simpler case of a zero or very small magnetic field 
in the case of strong thermal fluctuations was in fact 
comprehensively studied theoretically^ albeit in linear re- 
sponse only. In any superconductor there exists a critical 
region around the critical temperature \T — T c \ <C Gi-T c , 
in which the fluctuations are strong (the Ginzburg num- 
ber characterizing the strength of thermal fluctuations is 
just Gi~ 1CT 10 -I(r 7 for low T c , while Gi ~ lO^-lCT 1 
for HTSC materials). Outside the critical region and for 
small electric fields, the fluctuation conductivity was cal- 
culated by Aslamazov and Larkin^ by considering (nonin- 
teracting) Gaussian fluctuations within Bardeen-Cooper- 
Schrieffer (BCS) and within a more phenomenological 
Ginzburg-Landau (GL) approach. In the framework of 
the GL approach (restricted to the lowest Landau level 
approximation), Ullah and Dorseyi^ computed the Et- 
tingshausen coefficient by using the Hartree approxima- 
tion. This approach was extended to other transport 
phenomena like the Hall conductivity^ and the Nernst 
effect. 

The fluctuation conductivity within linear response 
can be applied to describe sufficiently weak electric fields, 
which do not perturb the fluctuations' spectrum™. Phys- 
ically at electric field, which is able to accelerate the 
paired electrons on a distance of the order of the co- 
herence length £ so that they change their energy by a 
value corresponding to the Cooper pair binding energy, 
the linear response is already inapplicable^. The result- 
ing additional field dependent depairing leads to devia- 
tion of the current-voltage characteristics from the Ohm's 
law. The non-Ohmic fluctuation conductivity was calcu- 
lated for a layered superconductor in an arbitrary elec- 
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trie field considering the fluctuations as noninteracting 
Gaussian one a 13 i 14 . The fluctuations' suppression effect 
of high electric fields in HTSC was investigated experi- 
mentally for the in-plane paraconductivity in zero mag- 
netic fieldi^— , and a good agreement with the theoreti- 
cal model o 13 i 14 was found. 

In this paper the nonlinear electric response of the 
moving vortex liquid in a layered superconductor under 
magnetic field perpendicular to the layers is studied us- 
ing the time dependent GL (TDGL) approach. The lay- 
ered structure is modeled via the Lawrence-Doniach dis- 
cretization in the magnetic field direction. In the moving 
vortex liquid the long range crystalline order is lost due 
to thermal fluctuations and the vortex matter becomes 
homogeneous on a scale above the average inter-vortex 
distances. Although sometimes motion tends to suppress 
the fluctuations, they are still a dominant factor in flux 
dynamics. The TDGL approach is an ideal tool to study 
a combined effect of the dissipative (overdamped) flux 
motion and thermal fluctuations conveniently modeled 
by the Langevin white noise. The interaction term in 
dynamics is treated in Gaussian approximation which is 
similar in structure to the Hartree-Fock one. 

Theoretically the nonlinear effects in HTSC have been 
addressed^. However the results of Ref. Il8j are different 
from our in this paper and we will sketch the difference 
below. 



is physically different 
18l believe that the two 



Firstly the model of Ref. |18, 
from ours. The authors in Ref. 
quantities, layer distance and thickness in the Lawrence- 
Doniach for HTSC are equal (apparently not the case 
in HTSC), while we consider them as two independent 
parameters. Another difference is we use so called self- 
consistent Gaussian approximation to treat the model 
while Ref. [TH used the Hartree-Fock approximation. 

A main contribution of our paper is an explicit form of 
the Green function incorporating all Landau levels. This 
allows to obtain explicit formulas without need to cutoff 
higher Landau levels. In Ref. Ila . a nontrivial matrix 
inversion (of infinite matrices) or cutting off the num- 
ber of Landau levels is required. Note that the exact 
analytical expression of Green function of the linearized 
TDGL equation in DC field can be even generalized also 
to AC field. The method is very general, and it allow us 
to study transport phenomena beyond linear response of 
type-II superconductor like the Nernst effect, Hall effect. 
The renormalization of the models is also different from 
Ref. [TH One of the main result of our work is that the 
conductivity formula is independent of UV cutoff (unlike 
in Ref. [Lsh as it should be as the standard |\l/| 4 theory 
is renormalizable. Furthermore Gaussian approximation 
used in this paper is consistent to leading order with per- 
turbation theory, see Ref. [l9| in which it is shown that 
this procedure preserved a correct the ultraviolet (UV) 
renormalization (is RG invariant). Without electric field 
the issue was comprehensively discussed in a textbook 
Kleinert 1 ^. One can use Hartree-Fock procedure only 
when UV issues are unimportant. We can also show, if 



there is no electric field, the result obtained using TDGL 
model and Gaussian approximation will lead the same 
thermodynamic equation using Gaussian approximation. 

The paper is organized as follows. The model is defined 
in Sec. II. The vortex liquid within the Gaussian approx- 
imation is described in Sec. III. The I-V curve and the 
comparison with experiment are described in Sec. IV, 
while Sec. V contains conclusions. 



II. THERMAL FLUCTUATIONS IN THE TIME 
DEPENDENT GL LAWRENCE-DONIACH 
MODEL 

To describe fluctuation of order parameter in lay- 
ered superconductors, one can start with the Lawrence- 
Doniach expression of the GL free energy of the 2D layers 
with a Josephson coupling between them: 
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where s' is the order parameter effective "thickness" and 
d! distance between layers labeled by n. The Lawrence- 
Doniach model approximates paired electrons DOS by 
homogeneous infinitely thin planes separated by distance 
d! . While discussing thermal fluctuations, we have to in- 
troduce a finite thickness, otherwise the fluctuations will 
not allow the condensate to exist (Mermin- Wagner the- 
orem). The thickness is of course smaller than the dis- 
tance between the layers (otherwise we would not have 
layers). The order parameter is assumed to be non- 
zero within s'. Effective Cooper pair mass in the ab 
plane is m* (disregarding for simplicity the anisotropy be- 
tween the crystallographic a and b axes), while along the 
c axis it is much larger-m c . For simplicity we assume 
a = aT™f(t -l),t= T/T™f, although this temperature 
dependence can be easily modified to better describe the 
experimental coherence length. The "mean field" critical 
temperature T™ 1 depends on UV cutoff, r c , of the "meso- 
scopic" or "phenomenological" GL description, specified 
later. This temperature is higher than measured critical 
temperature T c due to strong thermal fluctuations on the 
mesoscopic scale. 

The covariant derivatives are defined by D = V + 
i(27r/$o)A, where the vector potential describes con- 
stant and homogeneous magnetic field A = (—By, 0) and 
$0 = hc/e* is the flux quantum with e* = 2 |e|. The two 
scales, the coherence length £ 2 = ti 2 / (2m* aT c ) , and the 
penetration depth, A 2 = c 2 m*b' j '(4Tre* 2 aT c ) define the 
GL ratio K = A/£, which is very large for HTSC. In this 
case of strongly type-II superconductors the magnetiza- 
tion is by a factor k 2 smaller than the external field for 
magnetic field larger than the first critical field H c \ (T), 
so that we take B sa H. The electric current, J = J n +J s , 
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includes both the Ohmic normal part 



and the supercurrent 



2m' 



Jn = 0>,E, 



(2) 



(3) 



Since we are interested in a transport phenomenon, it is 
necessary to introduce a dynamics of the order param- 
eter. The simplest one is a gauge-invariant version of 
the "type A" relaxational dynamics 20 . In the presence of 
thermal fluctuations, which on the mesoscopic scale are 
represented by a complex white noise 2 ^, it reads: 



ay 

2m* 



1 8F, 



GL 



(4) 



where D T = d/dr— i(e* /h)$ is the covariant time deriva- 
tive, with <f> = — Ey being the scalar electric potential 
describing the driving force in a purely dissipative dy- 
namics. The electric field is therefore directed along the 
y axis and consequently the vortices are moving in the x 
direction. For magnetic fields that are not too low, we 
assume that the electric field is also homogeneous 2 ^. The 
inverse diffusion constant 7'/2, controlling the time scale 
of dynamical processes via dissipation, is real, although 
a small imaginary (Hall) part is also generally present 2 ^. 
The variance of the thermal noise, determining the tem- 
perature T is taken to be the Gaussian white noise: 



(C:(r,r)C m (r',r')> 



S 2 7 



7 TS{r - r')8(r - T')S nm . (5) 



Throughout most of the paper we use the coherence 
length £ as a unit of length and H C 2 = $o/27r£ 2 as a 
unit of the magnetic field. The dimensionless Boltzmann 
factor in these units is: 



n 



s 

It, 
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where the covariant derivatives in dimensionless units in 
Landau gauge are D x = J| - iby, D y = with b = 
B/H C 2, and the order parameter field was rescaled: "J 2 = 
(2aT™f /b')ip 2 . The dimensionless fluctuations' strength 
coefficient is 



u = \2Gm, 
where the Ginzburg number is defined by 

Gl = I (8eV£T c m Vc 2 fi 2 ) 2 



(7) 



(8) 



Note that here we use the standard definition of the 
Ginzburg number different from that in Ref. 23. The 



relation between parameters of the two models, the 
Lawrence-Doniach and the 3D anisotropic GL model, is 
d' = d^ c = d£,Hi s> = s £e = s£/7, where 7 2 = m c /m* 
is an anisotropy parameter. In analogy to the coher- 
ence length and the penetration depth, one can define a 
characteristic time scale. In the superconducting phase 
a typical "relaxation" time is tql = 7'£ 2 /2- It is con- 
venient to use the following unit of the electric field and 
the dimensionless field: Eql = H C 2^/ctgl,£ = E/Eql- 
The TDGL Eq. (j4j written in dimensionless units reads 

#V^^ + ^(2Vv ^ -VV l +l-VV l -l)- : ^^V^^ + |V^^lVn = Cn, 

(9) 

H = D T - l -D\ 

while the Gaussian white-noise correlation takes a form 

CT(r, t)U{v', t')\ = —5(r - v')6(r - r')8 nm . (10) 
/ s 

The covariant time derivative in dimensionless units 
is D T = -§- + ivby with v = £/b being the vortex 
velocity and the thermal noise was rescaled as = 
C n (2aT™f ) 3 / 2 /b' 1 / 2 . The dimensionless current density 
is J s = JglJs where 



(ii) 



with Jql — cH c2 / (2tt^k 2 ) being the unit of the current 
density. Consistently the conductivity will be given in 
units of (Jgl = Jgl/Egl = c 2 7' / (Att k 2 ) . This unit is 
close to the normal state conductivity a n in dirty limit 
superconductors^. In general there is a factor k of order 
one relating the two: a n — kaoL- 



III. VORTEX LIQUID WITHIN THE GAUSSIAN 
APPROXIMATION 

A. Gap equation 

Thermal fluctuations in vortex liquid frustrate the 
phase of the order parameter, so that (ip n (r,r)) = 0. 
Therefore the contributions to the expectation values 
of physical quantities like the electric current come ex- 
clusively from the correlations, the most important be- 
ing the quadratic one (^ n (r, t)^/>* (r', t')). In particular, 
(|?/>„(r, r)| 2 ) is the superfluid density. A simple approxi- 
mation which captures the most interesting fluctuations 
effects in the Gaussian approximation, in which the cu- 
bic term in the TDGL Eq. ©, \ip n \ 2 ip n , is replaced by a 
linear one 2 (IV'nl 2 ) V'n 



(12) 
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leading the "renormalized" value of the coefficient of the 
linear term: 



= - a/l + 2(|V„| 2 ) 



(13) 



where the constant is defined as = (1 — t — b)/2. The 
average { | V'n 1 2 ) is expressed via the parameter e below 
and will be determined self-consistently together with e. 
It differs slightly from a well known Hartree-Fock pro- 
cedure in which the coefficient of the linearized term is 
generally different (see^S for details). 

Due to the discrete translation invariance in the field 
direction z, it is convenient to work with the Fourier 
transform with respect to the layer index: 



ipn (r,r) 



2tt 



V^(r,r) - d^e^^r), 



(14) 



and similar transformation for £. In terms of Fourier 
components the TDGL Eq. (fT2"]) becomes 

H - ^ + ^[1 - cos{k z d)} + el V> fc , (r, r) = OT(r, t). 

The noise correlation is 

07* (r, r)OT(r', r')) = 47rcjt-d(T-v')S(T-r')S(k z -k' z ). 
I s 

. (16) 

The relaxational linearized TDGL equation with a 
Langevin noise, Eq. (|15[) . is solved using the retarded 
(G = for r < r') Green function (GF) G kz (r, r; r', t'): 



^n(r, r) 



The GF satisfies 



2 *l d dk z 



-ink z d 



dr' 



In 27T 

< ! dT'G kz (r,r;r',r')0~(r',T'). (17) 



H-^ + ^[l-co S (k z d)]+e\G, (r.r'.r- r') 



= <y(r-r')5(T-r'), 



(18) 



and is computed in the Appendix A. 

The thermal average of the superfluid density (density 
of Cooper pairs) is 

d t 2lxld dk 7 



hMr,r)|' 



2ut- 



2tt 



where 



/(£,t) 



x J dT'\G K (r-r',T-r'f 

uAb f°° f(e,r) 
27rs J T=T sinh(6r) ' 

exp ^! tanh (^l e -(^W). 
xe-^Io (2r/d 2 ) 



(19) 



(20) 



Here I (x) = (1/2%) e xcos6 d6 is the modified Bessel 
function. The first pair of multipliers in Eq. (|20[) is inde- 
pendent of the inter-plane distance d and exponentially 
decreases for r > (2e — b + v 2 ) , while the last pair of 
multipliers depends on the layered structure. The ex- 
pression (|T^|) is divergent at small r, so an UV cutoff t c 
is necessary for regularization. Substituting the expec- 
tation value into the "gap equation" , Eq. (|T3l) , the later 
takes a form 



-ah 



uAb f°° f(e,r) 
tts J T=Tc sinh(6r) 
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B. Renormalization 

In order to absorb the divergence into a "renormalized" 
value a r h of the coefficient , it is convenient to make an 
integration by parts in the last term for small r c : 



f(r) 
T—T c sinh(6r) 



ln[sinh(6r)] — 
dr 



/(g,T) 
cosh(frr) 



ln(6r c ). (22) 



Physically the renormalization corresponds to reduction 
of the critical temperature by the thermal fluctuations 
from to T c , The thermal fluctuations occur on the 
mesoscopic scale. The critical temperature T c is defined 
at e = 0, and v = 0, and at low magnetic field less than 
H c i = |p# In (k) (for a typical high T c superconductor, 
K ~ 50, H c i = 7.8 x 10~ 4 H c2 ), the superconductor is 
at Meissner phase, 6 = 0, leading to 



r c = Tr/|l + ^[ln(r c /d 2 )+ 7£ ]|, 



(23) 



where -fE — 0.577 is Euler constant, and Eq. (|2"Tj) can be 
rewritten as 



-a 



ln[sinh(6r)l — 



cosh(6r) 



+ — { 7£ -ln(M 2 )}, 



where a T h 
Gi = \( 



l-b-T/Tc 
2 



(24) 



t = T/T c , and u = \/2Gi% where 

z 2 K 2 ^T c j/c 2 h 2 ) 2 (T^ is now replaced by T c ). 
The formula is cutoff independent. In terms of energy UV 
cutoff A, introduced for example ir>ii, the cutoff "time" 
t c can be expressed as 



r c = l/(2e^A). 



(25) 



This is obtained by comparing a thermodynamic result 
for a physical quantity like superfluid density with the dy- 
namic result (see Appendix B). The temporary UV cutoff 
used is completely equivalent to the standard energy or 
momentum cutoff Lambda used in thermodynamics (in 
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which the time dependence does not appear). Physically 
one might think about momentum cutoff as more basic 
and this would be universal and independent of partic- 
ular time dependent realization of thermal fluctuations 
(TDGL with white noise in our case). Roughly (in phys- 
ical units) A ~ ef — h 2 k F /(2m*). In the next section 
we will discuss the estimate of T™^ using this value due 
to the following reason. For high-T c materials ordinary 
BCS is invalid and coherence length is of order of lattice 
spacing (the cutoff becomes microscopic) and therefore 
the energy cutoff is of order ep. Except the formula to 
calculate T" l f , all other formulas in this paper is inde- 
pendent of energy cutoff. 



IV. THE I-V CURVE 

A. Current density 

The supercurrent density, defined by Eq. (ITT1) . can be 
expressed via the Green's functions as: 



iu>t- 



2tt 



Gt 



r ,r 



r') 



x -^-G kz (r - t',t - t') + c.c. 
dy 

Performing the integrals, one obtains: 



4tts J t=q cosh(§) 2 ' 



3i 



(26) 



(27) 



where the function / was defined in Eq. (|20|) . Conse- 
quently the contribution to the conductivity is a s = jy/£- 
The conductivity expression (Eq. 27) is not divergent 
when expressed as a function of renormalized T c (the 
real transition temperature), so it is independent of the 
cutoff. This is considered in detail in section III. B and 
is indeed different from the Ref. 0. In physical units the 
current density reads 



Jy — O n E 



1 + 



ut_l f 00 f(e,r) 

4TTSkJ T=0 COSh(^l) 2 



(28) 



This is the main result of the present paper. We also 
considered the conductivity expression in 2D in linear 
response which do match the linear response conductivity 
expression derived in our previous workii. 



cut 



'2D 



2-1- 



2s 



Ansb \ \ b 
where ip is the polygamma function. 



(29) 



B. Comparison with experiment 

In this section we use physical units, while the dimen- 
sionless quantities are denoted with bars. The experi- 
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FIG. 1: Points are resistivity for different electric fields of 
an optimally doped YBCO in Ref. 0. The solid line is the 
theoretical value of resistivity calculated from Eq. (|30l) with 
fitting parameters (see text). The dashed line is the theo- 
retical value of resistivity in linear response with the same 
parameters. 



ment results of I. Puica et alZ, obtained from the re- 
sistivity and Hall effect measurements on an optimally 
doped YBa2Cu307_(5 (YBCO) films of thickness 50 nm 
and T c = 86.8 K. The distance between the bilayers used 
the calculation is d' = 11.68 A in Ref. |H The num- 
ber of bilayers is 50, large enough to be described by the 
Lawrence-Doniach model without taking care of bound- 
ary conditions. In order to compare the fluctuation con- 
ductivity with experimental data in HTSC, one cannot 
use the expression of relaxation time 7' in BCS which 
may be suitable for low-T c superconductor. Instead of 
this, we use the factor k as fitting parameter. The com- 
parison is presented in Fig. 1. The resistivity 

P=—^—, (30) 



er s + On 



(31) 



curves were fitted to Eq. (|30l) with the normal-state 
conductivity measured in Ref. Q to be a n — 1.9 x 10 4 
(ficm)^ 1 . The parameters we obtain from the fit are: 
H c2 (0) = T c dH c2 (T) /dT\ Tc = 190 T (corresponding to 
£ = 13.2 A), the Ginzburg-Landau parameter k — 43.6, 
the order parameter effective thickness s' = 8.5 A, and 
the factor k = <J n /&GL = 0.55, where we take 7 = 7.8 
for optimally doped YBCO in Ref. [26]. Using those pa- 
rameters, we obtain Gi = 9.32 x 10~ 4 (corresponding to 
u) = 0.136). The order parameter effective thickness s' 
can be taken to be equal to the layer distance (see in 
Ref. l27t ) of the superconducting CuO plane plus the co- 
herence length 2£ c = 2^ due to the proximity effect: 3.18 

A+2-yj- A= 6.9 A, roughly in agreement in magnitude 
with the fitting value of s' . 
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FIG. 2: The current-voltage curves calculated from Eq. (|28[) 
by using the parameters (see text) for different magnetic fields 
b = B/Hc2: 0.04 (1), O.f (2), 0.4 (3), 1.0 (4) at temperature 
t = 0.75. 




FIG. 3: The current-voltage curves calculated from Eq. (|28[) 
by using the parameters (see text) for different temperatures 
t = T/Tc: 0.2 (1), 0.3 (2), 0.4 (3), 1.0 (4) at magnetic fried 
6 = 0.5. 



We will now estimate T™^ for this sample. For the 
underdoped YBCO, the radius of the Fermi surface of 
YBCO was measured in Ref. EH k F = 0.7 A" 1 , while 
the effective mass is m* = 1.9m e . We will assume 
that the Fermi energy for underdoped YBCO of Ref. 
[H is sf = h 2 k F /(2m*) and is roughly the same for 
the optimal YBCO studied in this paper. The cut- 
off "time" in physical units is then, according to Eq. 
(j25]) . t c = 1.39 x KT 17 s. Equation (|23|) gives then 
T™* = 101.15 K. Using the parameters specified above 
we plot several theoretical FV curves. As expected the 
I-V curve shown in Figs. 2 and Figs. 3 has two lin- 
ear portions, the flux flow part for E <C Eql and the 
normal Ohmic part for E 3> Eql. In the crossover re- 
gion, E ~ Eql, a I-V curve becomes nonlinear due to 



destruction of superconductivity (the normal area inside 
the vortex cores increases to fill all the space). In Fig. 
2 the I-V curves are shown for different the magnetic 
fields, at a fixed temperature T = 0.75T C . At given elec- 
tric field, as the magnetic field increases, the supercur- 
rent decreases. When the magnetic field reaches H C 2, 
the I-V curve becomes linear. In Fig. 3 the I-V curves 
are shown for different temperatures, at a fixed magnetic 
field H = 0.5H C 2. At given electric field, as the tem- 
perature increases, the supercurrent decreases. When 
the temperature reaches T c , the I-V curve becomes lin- 
ear. With decreasing temperature the crossover becomes 
steeper. 



V. DISCUSSION AND CONCLUSION 

We quantitatively studied the transport in a layered 
type-II superconductor in magnetic field in the presence 
of strong thermal fluctuations on the mesoscopic scale 
beyond the linear response. While in the normal state 
the dissipation involves unpaired electrons, in the mixed 
phase it takes a form of the flux flow. Time depen- 
dent Ginzburg-Landau equations with thermal noise de- 
scribing the thermal fluctuations is used to describe the 
vortex-liquid regime and arbitrary flux flow velocities. 
We avoid assuming the lowest Landau level approxima- 
tion, so that the approach is valid for arbitrary values 
the magnetic field not too close to H c \ (T). 

Our main objective is to study layered high-T c mate- 
rials for which the Ginzburg number characterizing the 
strength of thermal fluctuations is exceptionally high, in 
the moving vortex matter the crystalline order is lost and 
it becomes homogeneous on a scale above the average 
inter-vortex distances. This ceases to be the case at very 
low temperature at which two additional factors make 
the calculation invalid. One is the validity of the GL 
approach (strictly speaking not far from T c (H)) and an- 
other is effect of quenched disorder. The later becomes 
insignificant at elevated temperature due to a very ef- 
fective thermal depinning. Although sometimes motion 
tends to suppress fluctuations, they are still a dominant 
factor in flux dynamics. The nonlinear term in dynam- 
ics is treated using the renormalized Gaussian approxi- 
mation. The renormalization of the critical temperature 
is calculated and is strong in layered high T c materials. 
The results were compared to the experimental data on 
HTSC. Our the resistivity results are in good qualitative 
and even quantitative agreement with experimental data 
on YBa2Cu307_5 in strong electric fields. 

Let us compare the present approach with a widely 
used Londons' approximation. Since we haven't ne- 
glected higher Landau levels, as very often is done in 
similar studie a 10 i 29 , our results should be applicable even 
for relatively small fields in which the London approxima- 
tion is valid and used. There is no contradiction since the 
two approximations have a very large overlap of applica- 
bility regions for strongly type-II superconductors. The 
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GL approach for the constant magnetic induction works 
for H >> H c i (T), while the Londons' approach works 
for H « H C 2 (T). Similar methods can be applied to 
other electric transport phenomena like the Hall conduc- 
tivity and thermal transport phenomena like the Nernst 
effect. The results, at least in 2D, can be in principle 
compared to numerical simulations of Langevin dynam- 
ics. Efforts in this direction are under way. 



The normalization is dictated by the delta function term 
in definition of the Green's function Eq. ([Tj 



Appendix B: Comparison with thermodynamics 



From TDGL, we obtain in the case v = 0: 



Acknowledgments 

We are grateful to I. Puica for providing details of ex- 
periments and to F. P. J. Lin, B. Ya. Shapiro for dis- 
cussions and encouragement, B. D. Tinh thanks for the 
hospitality of Peking University and D. Li thanks for 
the hospitality of National Chiao Tung University. This 
work is supported by NSC of R.O.C. #98-2112-M-009- 
014-MY3 and MOE ATU program, and D. Li is sup- 
ported by National Natural Science Foundation of China 
(grant number 90403002 and 10974001). 



|^n(r,r)|' 



2 \ ojtb f°° exp {- (2e ~ b + j?) t} I (%) 



2ns 



sinh(6r) 



(Bl) 

The supcrfluid density at b = and e = can be obtained 
by taking b and e to zero limit in the above equation: 



|^n(r,r)|' 



Lot 

27TS 7r=r c /d 2 



cx P {-2t}/ (2t) 



Performing the integration by parts, one obtains 



(B2) 



Appendix A: Derivation of the Green function of the 
linearized TDGL equation 



|^(r,r)| 2 ) ~ (r c /d 2 ) + lE } + (r c ) . (B3) 



In this appendix we outline the method for obtaining 
the Green function in strong electric field for the lin- 
earized equation of TDGL (fT5|) . The Green function is a 
Gaussian 



G K (r,rV 
where 



exp 



ib 

2 X (y + y) 



g k , (X,Y,t") 



(Al) 



X 2 +Y 2 
2/3 



-vX 



g K (X,Y, T ")=C kz (T")6(T")exp(- 

(A2) 

with X = x - x' - vt", Y = y-y',r" = t- t'. 6 (r") is 
the Heaviside step function, C and /3 are coefficients. 

Substituting the Ansatz (jAlj) into Eq. (fT8|) . one ob- 
tains following conditions condition: 

£-^ + y + ^[l-cos(M)] + i + %?=0, (A3) 



C 



1 



= 0. 



(A4) 



/3 2 j3 2 4 

The Eq. (|A4[) determines /3, subject to an initial condi 
tion /3(0) = 0, 



In the case without external electric field (or v = 
0), the equation obtained from TDGL shall ap- 
proach the thermodynamics result. In thermodynamics 
method, we shall evaluate the partition function Z = 
J Di}j n Di}* n e~ FGL / T where F GL /T is defined in Eq. ©. 
The superfluid density in the thermodynamic approach 
at the phase transition point 



\ipn(r,T)\ 



rfk 

{2irY s Jo 

2-K/d ^ 



dk z 



k2 l-cos(fc z d) 
2 + d 2 

~ ^{lnA + ln(2rf 2 )} + 0(A- 1 ), 

(B4) 

where A = k 2 max /2. 

The relation between the cutoff "time" t c and energy 
UV cutoff A is obtained by comparing Eq. (|B3|) with Eq. 
(IB41) 



f3 = - tanh (br"/2) . 



while Eq. (|A3|) determines C: 



C = — exp 

47T 



x sinh 



br" 



(A5) 



S-5 + Y + ^[l-CO S (M)]) r" 



(A6) 



1 



2AeT E ' 



(B5) 



We also remark that in thermodynamic approach, if 
we use the Gaussian approximation, we will get the exact 
same equation derived in Eq. (1241) without electric field 
derived from TDGL after using Eq. (|B5[) . 
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